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Abstract 

We investigate in this paper the cosmological evolution of a dark energy model with two 
scalar fields where one of the scalar has canonical kinetic energy and another scalar has 
negative kinetic energy term. For such a system with exponential potentials we find that 
during the evolution of the universe the equation of state w changes from w > — 1 to 
w < — 1, which is consistent with the recent observations. A phase-plane analysis shows 
that the "phantom" -dominated scaling solution is the stable late-time attractor of this type 
of models. 
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Scalar fields play an important role in modern cosmology. The dark energy can be 
attributed to the dynamics of a scalar or multi-scalar fields, for instance the quintessence 
PI 12], which convincingly realize the present-day cosmic acceleration by using late-time 
attractor solutions, in which the scalar fields mimic the perfect fluid in a wide range of 
parameters. For the detail studies on the models much attention has been drawn to the 
case of exponential potentials. The exponential potentials allow the possible existence of 
scaling solutions in which the scalar field energy density tracks that of the perfect fluid (so 
that at late times neither field is negligible). In particular, a phase-plane analysis of the 
spatially flat FRW models showed that these solutions are the unique late-time attractors 
whenever they exist jSl E]. Moreover, exponential potentials often appear naturally in 
models after compactification in string/M-theory. 

By fitting the recent SNe la data, marginal (2a) evidence for w(z) < — 1 at z < 0.2 has 
been found |5\. To obtain w < — 1, phantom field with a negative kinetic term may be a 
simplest implementing, in which the weak energy condition is violated, and can be regarded 
as one of interesting possibilities describing dark energy |H]. The physical background 
for phantom type of matter with strongly negative pressure may be looked for in string 
theory [2]. Phantom field may also arise from a bulk viscous stress due to the particle 
production [H] or in higher-order theories of gravity ]H], Brans-Dicke and non- minimally 
coupled scalar field theories ^0] • The cosmological models which allow for phantom matter 
appear naturally in the mirage cosmology of the braneworld scenario ^T] and in k-essence 
models |12j . In spite of the fact that the field theory of phantom fields encounters the 
problem of stability which one could try to bypass by assuming them to be effective fields [131 
ITlj . it is nevertheless interesting to study their cosmological implication. Recently, there are 
many relevant studies on phantom energy ^H] and the primordial perturbation spectrum 
from various phantom inflation models [16] . 

Furthermore, the analysis on the properties of dark energy from the recent observations 
mildly favor models with w crossing —1 in the near past. However, neither quintessence nor 
phantom can fulfill this transition. But in an universe with a quintessence and a phantom 
this case can be realized easily. This implement of dark energy, called as quintom, has 
been proposed in Ref. [Hj, and in some cases providing a better fit to the data than the 
more familiar models with w > — 1. The quintom might be able to preserve the tracking 
behavior of the quintessence and needs less fine-tuned in the early universe compared with 
the phantom field, furthermore predict a interesting feature in the evolution and fate of the 
universe ^Hj- In this paper we study the the quintom model above with an exponential 
potential and investigate the existence and stability of cosmological scaling solutions in 
the context of spatially homogeneous cosmological models. Our phase-plane analysis of 
the spatially flat FRW models shows that the phantom-dominated scaling solution is the 
unique late-time attractor and there exists a transition from w > —1 to w < — 1, then to 
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a constant related to the slope of the phantom field potential at late times. We will also 
discuss the physical consequences of these results. 

We consider a toy model which contains a negative-kinetic scalar field and a normal 
scalar field a with generic exponential potential, described by the action: 

3 = J d 4 xV=g - \<T%WA + \g^d,od u a + V(<j>, a) + £ m ) , (1) 

where k 2 = 8tiGn is the gravitational coupling and Cm represents the Lagrangian density of 
matter fields. The homogeneous fields and a in a spatially flat FRW cosmological model 
can be described by a fluid with an effective energy density p and an effective pressure P 
given by 

P = -\tf + \a 2 + V^a) } (2) 
P = -^ 2 + ^ 2 -V(0,a). (3) 
The corresponding equation of state parameter is now given by 

-0 2 + a 2 -2\/(0,tx) 

W = : . (4) 

-(p 2 + & 2 + 2V((f>, a) 

For a model with a normal scalar field, the equation of state w > — 1. The toy model of 
a phantom energy component with a negatice kinetic term possesses an equation of state 
w < —1. In our model, Eq.(|4"J) implies w > — 1 when & > and w < — 1 when a < <fi. We 
assume that there is no direct coupling between the phantom field and the normal scalar 
field with such a potential 

V{cj>, a) = VM) + V a (a) = V^ e~ X ^ + V a0 e"^, (5) 

where and X a are two dimensionless constants characterising the slope of the potential 
for (p and a respectively. Further we assume > and X a > since we can make them 
positive through — > —0 and a — > —a if some of them are negative, respectively. The 
evolution equations of the fields and the fluid for a spatially flat FRW model with Hubble 
parameter H is 

+ 3iJ0 = 0, (6) 

cr + 3iJcr H = 0, (7) 

da 

p 7 + 3#(p 7 + P 7 ) = 0, (8) 

where p 7 is the density of fluid with a barotropic equation of state P 7 = (7 — l)p 7 , where 
7 is a constant, < 7 < 2, such as radiation (7 = 4/3) or dust (7 = 1). The Fridemann 
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constraint equation is 

H 2 = y (-^ 2 + VM + \° 2 + V a (a) + p 7 ) • (9) 



Defining five dimensionless variables 



1V "7& ' s - = 7S^' (10) 
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the evolution equations (JEJ)-© can be rewritten as an autonomous system: 
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= -3x [1 + - x a - -z j - A — , (11) 

v<t> = %y<t> [^ x l + + \ z2 - A«^^l , (12) 

x' a = -3x a (l + xl-xl- ^-z 2 ) +X a -^-yl, (13) 
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2/a = 3yJ-xl + xl + ^z 2 -\ a ^§-x a ), (14) 



= 3^-xJ + 4 + p-|), (15) 

where a prime denotes a derivative with respect to the logarithm of the scale factor, N = 
In a, and the Fridemann constraint equation Q becomes 

-xl + yl + xl + yl + z 2 = l. (16) 

We will restrict our discussion of the existence and stability of critical points to expanding 
universes with H > 0. Critical points correspond to fixed points where x'^ = 0, y'^ = 0, 
x' a = 0, y' a = and z 1 = 0, and there are self-similar solutions with 

= 3x1 - 3x l - Tp- ( 17 ) 

This corresponds to an expanding universe with a scale factor a(t) given by a oc t p , where 

2 

^ —6a;? + 6x 2 + 3^z 2 ^ 
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Table 1: The properties of the critical points in a spatially flat FRW universe containing a 
phantom field and a normal scalar field with exponential potentials. 



The system (IllJl - ffTo'j) has at most one two-dimensional hyperbola K embedded in five- 
dimensional phase-space corresponding to kinetic-dominated solutions, a fixed point P 
which is a phantom-dominated solution, a fixed point S which is a scalar-dominated solu- 
tion, a fixed point F which is a fluid-dominated solution, and a fixed point T which is a 
fluid-scalar-dominated solution listed in Table 1. 

In order to study the stability of the critical points, using the Friedmann constraint 
equation (fTtjj) we first reduce Eqs. (fTT|) - (fTo]) to four independent equations. Substituting 
linear perturbations x^ — > x^ + Sx^,, y^ — + y^ + Sy^, x a — > x a + 5x a and y a — > y a + 5y a about 
the critical points into the four independent equations, to first-order in the perturbations, 
gives the evolution equations of the linear perturbations, which yield four eigenvalues mj. 
Stability requires the real part of all eigenvalues to be negative. 

K: These kinetic-dominated solutions always exist for any form of the potential, which 
are equivalent to stiff-fluid dominated evolution with a oc t 1 / 3 irrespective of the nature 
of the potential. The linearization of system (fTTj) -()16J ) about these fixed points yields four 
eigenvalues 

mi = 0, TO2 = 3, 777-3 = 3(2 — 7), 7774 = 3(1 ± — p=), 

V6 

where we use upper /lower signs to denote the two distinct cases of x a = T\J^- + Thus 
the kinetic-dominated solutions are always unstable. 

P: The phantom-dominated solution exist for any A^, and \ a . The power-law exponent, 
p = — 2/A|, depends on the slope of the potential V</>. The equation of state becomes 
w = —1 — A|/3. The linearization of system (fTTf - (PT)|) about this fixed point yields four 
eigenvalues 

A 2 1 

mi = ^, 777 2 = 777 3 = --(6 + A|), 777 4 = - (37 + A|) , 

which indicate that the solution is stable. 



4 



S: The scalar-dominated solution exist for A 2 < 6. The power-law exponent, p = 2/\„, 
depends on the slope of the potential V a . The equation of state becomes w = — l + A^/3. 
The linearization of system (|11 )) -([16 )) about this critical point yields four eigenvalues 

A 2 1 
mi = y, rn 2 = m 3 = -(A 2 - 6), m 4 = A 2 - 37, 

which indicate that the solution is unstable. 

F: The fluid-dominated solution exists for any form of the potential, corresponding to a 
power-law solution with p = 2/3'j. The linearization of system (|11 )) -([16 )) about this critical 
point yields four eigenvalues 

37 37 
m 1 = m 2 = — , m 3 = m 4 = — - 3, 

which indicate that the solution is unstable. 

T: The scalar-fluid-dominated solution exist for a potential with A 2 > 37. The power- 
law exponent, p = 2/37, is identical to that of the fluid-dominated solution, depends only 
on the barotropic index 7 and is independent of the slope A CT of the potential V a . The 
linearization of system (JTTJ)- (JTOJ) about the fixed point yields four eigenvalues 

mi 
m 2 



m 3 



m 4 

which indicate that the solution is unstable. 

The case with the evolution of the state equation parameter w crossing —1, a scenario of 
quintom has be shown to be favored mildly by the recent observations. In this paper we have 
discussed a possible realization of it, in which both quintessence field and phantom field 
are introduced. We have presented a phase-space analysis of the evolution for a spatially 
flat FRW universe containing a barotropic fluid and phantom-scalar fields with exponential 
potentials and shown that the phantom-dominated scaling solution is the stable late-time 
attractor. The energy density of the phantom field dominates at the late time. The reason 
for this behavior is that the energy density of the phantom field increases while those of the 
normal scalar field and the barotropic fluid decreases as the universe evolves. Our numerical 
studies indicate that the state equation parameter w changes from above —1 to below — 1 
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w 




and tends to be —(1 + A|/3) as shown in Figure I. Moreover, for various selections of 
potentials, we find that the state equation parameter w changes from above —1 to below 
— 1 and tends to be —1 as shown in Figure 2, and from below —1 to above —1 and tends to 
be —1 in Figure 3. we have assumed that there is no direct coupling between the phantom 
field and the normal scalar field in this paper. 



Acknowledgements 

This project was in part supported by National Basic Research Program of China under 
Grant No.2003CB716300 and also by NNSFC under Grant No.10175070. 



References 

[1] B.Ratra and P.J. E.Peebles, Phys.Rev. D37 (1988) 3406; C.Wetterich, Nucl.Phys. 
B302 (1988) 302. 

[2] I.Zlatev, L.Wang and P. J.Steinhardt, Phys.Rev.Lett, 82 (1999) 896; P.J.Steinhardt, 
L.Wang and I.Zlatev, Phys.Rev. D59 (1999) 123504. 



6 




w 



-0 . 92 








-0 . 94 








-0 .96 








-0 . 98 








-1 . 02 


J r 4 


6 


8 10 


-1 . 04 









t 



Figure 3: The evolution of the effective equation of state of the phantom and normal scalar 
fields with V(<f>, a) = e" A ^ 2 + V a0 e~ x ^ a \ 
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